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We investigate the effects of charge independence and charge symmetry breaking in neutron- 
rich matter. We consider neutron and proton properties in isospin-asymmetric matter at normal 
densities as well as the high-density neutron matter equation of state and the bulk properties of 
neutron stars. We find charge symmetry and charge independence breaking effects to be very small. 



I. INTRODUCTION 

Charge independence (CI) of the nuclear force signi- 
fies that, in a state of total isospin equal to one (the 
only accessible to two neutrons or two protons), the 
neutron- proton (np), neutron- neutron (nn), and proton- 
proton (pp) interactions are equal, after electromagnetic 
effects have been removed. Charge independence break- 
ing (CIB) refers to slight violations of this property. On 
the other hand, charge symmetry breaking (CSB) refers 
only to differences between the nn and the pp interac- 
tions. Clearly, CSB also breaks CI, while the opposite is 
not necessarily true. Basically, charge independence is a 
statement of invariance under rotations in isospin space, 
whereas CS is a special case of CI and amounts to in- 
variance under a rotation from T z — 1 to T z = — 1. A 
thorough review of charge dependence can be found in 
Ref. pp. 

At a fundamental level, charge dependence of the nu- 
clear force is understood in terms of differences between 
the masses of up and down quarks and electromagnetic 
interactions among the quarks. At the hadronic level, 
CSB manifests itself through mass differences between 
hadrons of the same isospin multiplet, meson mixing, and 
other mechanisms such as meson-photon exchanges. 

Both CSB and CIB are seen most clearly in the 1 So 
scattering lengths, which have the values [THE] 

a nn = -18.9 ± 0.4 fm , a pp = -17.3 ± 0.4 fm , 

a np = -23.74 ± 0.02 fm . (1) 

These values for the scattering lengths imply that the nn 
interaction is slightly more attractive than the pp, and 
the (T=l) np interaction is more attractive than both pp 
and nn. The difference between a nn and a pp is a measure 
of CSB, whereas 

A 2 iP'nn ~t~ Qpp) Qnp (2) 

is typically taken as a measure of CIB. Thus, the sin- 
glet scattering length of nucleon-nucleon (NN) scattering 
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shows small CSB (about 10%) and larger CIB. It has been 
shown that both nucleon mass splitting [7] and p° — us 
mixing [5] can fully explain charge symmetry breaking in 
the singlet scattering length. Thus, different observables 
need to be identified which might reveal sensitivity to the 
mechanism used to describe CSB. 

In this paper, CSB/CIB are introduced in two differ- 
ent ways. In one case, we use the high-precision, charge- 
dependent NN potential known as the "CD-Bonn" po- 
tential [5] • The latter is based upon the Bonn full model 
[lOj and reproduces the predictions of the comprehensive 
model for CSB and CIB in all partial waves with J < 4. 
The other model uses p — ui mixing for CSB, with the 
mixing Lagrangian as in Ref. [TTJ. 

Within our recent efforts to investigate diverse as- 
pects of the nuclear equation of state (EoS) [12], our 
main purpose in this paper is to explore the effects of 
charge dependence and asymmetry of the nuclear force 
in a highly neutron-rich environment, such as isospin- 
asymmetric nuclear matter (IANM), and how the effect 
evolves as a function of density and/or isospin asymme- 
try. 

Isospin-asymmetric nuclear matter simulates the inte- 
rior of a heavy nucleus with unequal densities of protons 
and neutrons. The equation of state of (cold) IANM is 
then a function of density as well as the relative con- 
centrations of protons and neutrons. The recent and 
fast-growing interest in IANM stems from its close con- 
nection to the physics of neutron-rich nuclei, or, more 
generally, asymmetric nuclei, including the very "exotic" 
ones known as "halo" nuclei. The equation of state of 
IANM is also the crucial input for the structure equa- 
tions of compact stars, and thus establishes the connec- 
tion between nuclear physics and compact astrophysical 
systems. Hence, in-depth studies of IANM are important 
as well as timely. 

In symmetric nuclear matter, CSB would cause a split- 
ting of the single-particle potential, AU = U p — U n . 
Some calculations within the mean field approximation 
|13j find that CS is gradually restored in nuclear mat- 
ter in /^-equilibrium, whereas CSB becomes strongly en- 
hanced away from /3-equilibrium because of the large nu- 
cleon isovector density. Medium effects on CSB in neu- 
tron matter were investigated in neutron matter |14j and 
found to reduce the mass of a neutron star (composed of 
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neutrons only) by as much as 35% as compared to the one 
resulting from a standard relativistic calculation with no 
CSB effects. 

In the next Section, we will describe how CSB and CIB 
effects are included in the IANM environment. We will 
then present and discuss some results, with a particu- 
lar eye on selected isovector properties, specifically the 
single-proton and single-neutron potentials, the neutron- 
proton mass splitting, and the symmetry potential. Neu- 
tron star masses and radii will also be shown. Conclu- 
sions and future plans are summarized in the last Section. 



II. DESCRIPTION OF THE CALCULATION 

A. The CSB and CIB mechanisms through mass 
splitting 

CD-Bonn is a high-precision charge-dependent one- 
boson-exchange potential whose charge dependence is 
based upon the predictions of the Bonn full model [10] . 
The charge symmetry breaking due to nucleon mass split- 
ting was investigated in Ref. [7] with the Bonn full model. 
Considerable CSB was found to be generated by 2w and 
7r — p diagrams. In addition to the well-known CSB in 
the singlet scattering length, non-negligible CSB effects 
were observed in higher partial waves, especially P and 
D waves. Apparently, those are crucial for a quantitative 
account of the Nolen-Shiffer [15 anomaly seen in the en- 
ergies of neighboring mirror nuclei |16j . 

Concerning CIB, pion mass splitting was found to be 
a major cause for it, in the one-pion as well as multi- 
pion exchanges. In fact, the contribution from multi-pion 
exchanges to CIB is about 50% of the one originating 
from one-pion exchange [T71 [TB] . 

CD-Bonn consists of three different potentials, V pp , 
V nn , and V pn , for the pp, nn, and pn cases, respectively. 
First, V pp is constructed through hts to the pp data. 
Next, starting from V pp and replacing the mass of the 
proton with the mass of the neutron, the nn singlet scat- 
tering length is reproduced. The V np in the T=l channel 
is obtained using the proper charge-dependent one-pion 
exchange plus a slight adjustment of the aNN coupling 
to account for the remaining CIB. The predictions by 
CD-Bonn are 



Here we use the Brueckner-Hartree-Fock (BHF) de- 
scription of nuclear matter. The scattering equation for 
two nucleons in nuclear matter reads, schematically 



G, 



(4) 



where G\j is the in-medium reaction matrix (ij=nn, pp, 
or np). Qij is the Pauli operator, which prevents scat- 
tering to occupied nn, pp, or np states. In the case of 
CD-Bonn, the scattering equation (in free space) is de- 
rived from the Blankenbecler and Sugar |19j reduction of 
the of the four-dimensional Bethe-Salpeter equation [2D] , 
and is formally identical to the non-relativistic Lippman- 
Schwinger equation. 

The goal is to determine self-consistently the nuclear 
matter single-particle potential which, in IANM, will be 
different for neutrons and protons. We have, for neu- 
trons, 



U n = U, 



and for protons 



np 



Up — Up n 



u, 



(5) 



(6) 



48.9680 im,a p 



47.4602 fm 



where each of the four pieces on the right-hand-side of 
Eqs. (5 6) signifies an integral of the appropriate G- 
matrix (nn, pp, or np/pn). The splitting of the single- 
particle potential manifest in Eqs. (5-6) is due to the pres- 
ence of two different Fermi momenta for neutrons and 
protons, regardless charge effects. In addition, charge 
symmetry and charge independence breaking can then 
be incorporated in a natural way using the appropriate 
potential in each of the integrals above, whereas in stan- 
dard charge-independent calculations the same potential 
(namely, the np potential) is used for all nucleon pairs. 
Thus, for an arbitrary level of isospin asymmetry, the ef- 
fects we include are both CIB and CSB effects. On the 
other hand, in nearly pure neutron matter, where the 
proton Fermi momentum approaches zero, U np in Eq. (5) 
and U pp in Eq. (6) will approach zero. The U pn term in 
Eq. (6) becomes the (very attractive) potential of a pro- 
ton impurity, but it has essentially zero weight in the final 
averaging of neutrons and protons due to the vanishing 
proton Fermi momentum. Therefore, using the CD-Bonn 
nn potential in the second term on the right-hand side of 
Eq. (5) instead of the np potential (as usually done) is to 
say that the mechanism under consideration is primarily 
CIB. 



a np = -23.7380 fm . (3) 



B. Application to isospin-asymmetric nuclear 
matter 

The framework we use to calculate the EoS of IANM 
P^] naturally lends itself to the application of potentials 
which are different for different nucleon pairs. 



III. RESULTS AND DISCUSSION 
A. Some neutron and proton properties in IANM 

All charge-dependent effects shown in this subsection 
are included through the use of the appropriate version 
of the CD-Bonn potential in Eqs. (5-6). 

We begin with the splitting of the neutron and proton 
effective masses in IANM, see Fig. I, as a function of the 
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FIG. 1: Neutron (green and blue lines) and proton (red and 
pink lines) effective masses in IANM as a function of the neu- 
tron excess parameter and for fixed density. The Fermi mo- 
mentum is equal to 1.4 fm _1 and 2 fm _1 in the upper and 
lower frame, respectively. The green dashed and the pink 
dotted lines show the predictions with CSB/CIB. 

neutron excess parameter, a = p " , Pp , in terms of neu- 

PnrPp 

tron and proton densities. The large splitting between 
the effective masses, which vary almost linearly and in 
opposite directions, is of course due to the presence of (in- 
creasing) isospin asymmetry in the medium. The much 
smaller effect, noticeable only on the right side mostly 
through the differences between the green and the blue 
lines, is due to charge dependence. Comparison between 
the predictions on the left side and on the right side of 
the figure suggests that these effects increase slightly with 
total density, and that they are larger for neutrons. 

Figures 2 and 3 show charge-dependent effects on the 
single-neutron potential as a function of the momentum 
and fixed density. Comparison between the predictions 
on the left side and on the right side of both figures in- 
dicates that charge-dependence effects increase only very 
little as the neutron excess parameter changes from 0.5 
to nearly 1. Taking CIB properly into account makes U n 
more repulsive because the nn interaction is less attrac- 
tive than the np one. 
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FIG. 2: The single-neutron potential in IANM as a function of 
momentum (in units of the Fermi momentum) for fixed den- 
sity and neutron excess parameter. In both frames the Fermi 
momentum is equal to 1.4 fm _1 , whereas the neutron excess 
parameter is equal to 0.5 in the upper frame and nearly 1 (al- 
most all neutrons) in the lower frame. The pink dotted line, 
almost indistinguishable from the solid red in the left panel, 
corresponds to the calculation including charge-dependent ef- 
fects. 



The symmetry potential, defined as 



plays a crucial role in separating the dynamics of neu- 
trons and protons in collisions of asymmetric systems. It 
can be related to the isovector component of the nuclear 
optical potential, and thus is a valuable quantity when 
seeking constraints to the EoS of IANM. In Fig. 4, we 
show the symmetry potential as a function of momentum 
and for two different values of the isospin asymmetry. Be- 
cause U n and U p are approximately linear with respect 
to a, the symmetry potential should be nearly indepen- 
dent of a, which is confirmed in Fig. 4. Again, the dotted 
line displays predictions obtained with charge-dependent 
effects. We observe that, although small, the difference 
is noticeable, particularly at the lowest momenta. 
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FIG. 3: As in Fig. 2, but with a Fermi momentum equal to 2 
fm" 1 . 

B. Masses and radii of neutron stars 

From the results shown in the previous Section, one 
might conclude that CSB and CIB effects in nuclear mat- 
ter, even if very neutron rich, are quite small. Neverthe- 
less, it is insightful to take the (charge-dependent) EoS to 
the very high densities probed by compact stars and cal- 
culate neutron star bulk properties. We will show these 
predictions as obtained with both the mass splitting and 
the p — u) mixing models. 

In the latter, the charge asymmetric nuclear force is 
evaluated by adding a "mixed" p — u exchange. The 
mixing matrix element was extracted in Ref. [5] from the 
e + e~ — > 7r + 7r~ annihilation process and found to be 

< p°\H m \uj >= -0.00452 ± 0.0006 GeV 2 . (8) 

Within the p — uj mixing model, the CSB potential is de- 
veloped starting with the pp CD-Bonn potential, which is 
the most accurate due to the large base of high-quality pp 
data, and then constructing the nn potential by adding 
to the pp potential a term that contains the difference 
Vpu = V pu (nn) - V pul {pp) generated by p - u> mixing. 
The reader is referred to Ref. [TT] for more details. 



FIG. 4: The symmetry potential without (solid red) and with 
(dotted pink) CSB/CIB effects as a function of momentum 
and fixed Fermi momentum, equal to 2 fm -1 . The neutron 
excess parameter is equal to 0.5 in the upper frame and nearly 
1 in the lower one. 



All predictions shown in this Section are calculated 
with: 1) No consideration of charge dependence or charge 
asymmetry (curve labeled "No CD" in Figs. 5-8); 2) CIB 
and CSB effects as included in CD-Bonn (the curve la- 
beled "CD, mass splitting" in Figs. 5-8); and 3) p — oj 
mixing for the CSB sector. At this point, it is useful to 
recall the comments made following Eqs. (5-6) concern- 
ing which kind of charge effects we may potentially detect 
in a systems of only neutrons. 

In Figs. 5-6 we see CIB effects on the energy per par- 
ticle and the pressure in neutron matter. The effect of 
CIB is small but noticeable, particularly at the higher 
densities, However, the predictions from the two models 
are essentially indistinguishable. The inclusion of CIB 
renders the EoS slightly more repulsive. This is reason- 
able, as the interaction between two identical nucleons 
described in terms of a charge-independent NN potential 
(typically np), is less repulsive than one designed specif- 
ically for two neutrons. 

The impact of CIB on the mass-radius relation of a 
neutron star (comprised of only neutrons) is shown in 
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Fig. 7. We note in passing that, in spite of the attractive 
nature of CD-Bonn, and, particularly, the use of the con- 
ventional BHF calculation of the EoS, one might have ex- 
pected even smaller maximum masses than those shown 
in Fig. 7. However, we must keep in mind that only neu- 
trons are involved here (i.e. the T=l potential). The 
over-attraction expected from CD-Bonn comes mainly 
from the T=0 potential. 

Back to CIB considerations, the maximum mass 
changes from 1.58 to 1.6 solar masses due to CIB, a varia- 
tion of only 1.3%. In Fig. 8, the mass of the star is shown 
as a function of the central energy density for the three 
models, showing that the central density corresponding 
to the maximum mass is approximately the same in all 
cases. 

In closing this Section, we take note of the findings 
from Ref. [T3], where a large reduction of the star mass, 
as much as 35% around normal density, was reported 
as a consequence of CSB. In that work, p — to mixing 
takes place through a baryon loop, as initially proposed 
by Piekarewicz ans Williams [31], and is subjected to 
medium effects. In contrast, we attribute p — uj mixing 
to electromagnetic coupling, see Eq. (8). 

IV. SUMMARY AND CONCLUSIONS 



FIG. 6: The pressure in neutron matter for the three models 
discussed in the text. 
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FIG. 7: The mass of a neutron star vs. the radius for the 
three models under consideration. 



Symmetries and the mechanisms responsible for their 
breaking are issues of fundamental importance as they 
help our understanding of the nuclear force at its very 
core. 

In this paper, we have investigated the impact of CIB 
and CSB of the NN interaction on the splitting between 
neutron and proton properties in IANM. We find such 
impact to be very small, even at very high neutron rela- 
tive concentrations. Given the chief role of the symmetry 
potential in collisions of isospin-asymmetric systems, it 
may be useful to keep in mind that there is some small 
but noticeable sensitivity at low momentum. 

We then moved to high-density neutron matter and 
calculated neutron star masses and radii. Consistently 
with the EoS becoming slightly more repulsive with con- 



siderations of CIB, the star mass is found to increase by 
a few percent. 

We noted that, from calculations where (density- 
dependent) nucleon-antinucleon loops allow the transi- 
tion between isoscalar and isovector mesons, a strong en- 
hancement of CSB is predicted in neutron matter. This 
is definitely an important issue to be further explored, 
as a large CSB effect in stellar matter would impact the 
cooling mechanism of the star. 

In the near future, we plan to explore if and how the 
presence of relativistic effects impact CIB and CSB in 
stellar matter. To that end, we will produce a charge- 
dependent relativistic potential (such as a "charge- 
dependent Bonn B" ) suitable for applications within the 
relativistic Dirac-Brueckner-Hartree-Fock approach. 
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FIG. 8: The mass of a neutron star vs. the central density 
for the three models under consideration. 
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